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APstract 

The methods used to analyte the aerodynamic 
performance of V/SIOL Inlets at the NaSa lewis Ke- 
search Center is briefly described. Kecent exten- 
sions and applications of the method are empha- 
slted. They Include the specification of the Kutta 
condition tor a slotted inlet, the calculation of 
suction and tangential blowing for boundary layer 
control, and the analysis of auxiliary inlet geome- 
tries at angles of attack. A comparison Is made 
with experiment tor the slotted inlet. Finally, an 
optimum diffuser velocity distribution is developed. 

Nomenc lature 

Cf skin friction coefficient 

U tan d I aw ter 

I inlet length 111. 3048 m) 

M Nach number 

lii mass flow rate 

5 Surface distance 

V veto. . . 

X K-coordinate 

"x nond linens lona I distance from the 

start of diffusion norm.ili;ed 
by the total diffusion length 
a angle of attack 

a inlet yaw angle 

e circumferential angle (e • 0* at wind- 

ward plane) 

^iitm^rj_g^ls: 

c control station 

de diffuser exit 

e edg«* of the boundary layer 

inlet inlet 

j billwing jet 

m,ix nuiximum 

ref referetue 

suction suction 

T throat 

» free stream 

Introduct ion 

In recent years, many different airfranx'/xn- 
gine configurations have been proposed for V/blUl 
aircraft, bonn' of the proposed conf igurat ions im- 
posed rather severe flow ccmditions on the propul- 
sion system inlet, for example, the approach and 
takeoff flight paths of a tilt nacelle V/blOl air- 
craft may result in inlet angles of attack up to 
UO*. A ncijor concern for the designer at these 
conditions is the possibi 1 ity of inlet internal 
flow separation, beparation tree flow is desired 
to minimi/e thrust loss, minimi/e tan blade stress, 
and prevent engine stall. Keliable theoretical 
methods of inlet flow analysis are desired to in- 
terpret and augment the results of wind tunnel 
testing, ihe wthods should be able to calculate 
the potential and boundary layer flows in inlets of 
arbitrary geowtry and flow conditions. 


buch methods ot analysis have been developed 
over the past several years at NASA lewis Research 
Center. Ihey consist of a series of computer pro- 
grams documented In Kefs. 1 to b. Comparisons with 
experimental results are presented In Kefs. 7 to 
U. Since these reports, the programs have been 
extended and applied to more difficult inlet prob- 
lems. 

The present paper will briefly describe the 
basic method of analysis. The major emphasis, how- 
ever, will be on presenting the recent extensions 
and applications. The topics covered In this paper 
are: the flow about a slotted Inlet; the flow In 

an auxiliary Inlet; the analysis of suction and 
blowing boundary layer control; and the developront 
of an optimum diffuser velocity distribution. 

Bas ic Method o f A nalysis 

Ihe basic problem to be solved 1$ to calculate 
the compressible viscous flow in inlets of arbi- 
trary geometry and operating conditions. A series 
of ccvnputer programs devt*lo;Hfd at NAbA lewis Re- 
search Center are used to solve this problem. A 
flow chart depicting tne sequence for using these 
programs is presented in Fig. 1 with the oasic pro- 
gr.ims on the left, and recent extensions on the 
right. All programs start with the geometry pro- 
gr.im, upper left-hand block, which creates the dis- 
crete control points for each geoim*tr1c configura- 
tion. Ihen the incompressible potential flow pro- 
gram is used to calculate the basic solutions to 
the problem. Ihese basic solutions are combined 
into a solution that satisfies the Inlet operating 
conditions of freestream velocity, angle of attack, 
•ind inlet mass flow. Next, the incompressible flow 
IS corrected for ccimpressihi 1 1 ty effects. Ihe com- 
pressible potential flow solution is then used as 
an input to the boundary layer program which calcu- 
lates the laminar, transltum and turbulent boundary 
layer characteristics, and predicts flow separation. 

Iwo Iteration loops are available as shciwn to 
the left In Fig. 1. Ihe first adds the displace- 
wnt thicknesc lo the gc'ometry to improve the ac- 
curacy of the pr.tential flow and boundary layer 
calculat Ions. Ihe second incorporates an automatic 
angle ot attack sweep to find the separation bound- 
ary of an inirt in one uninterrupted computer run. 

Ihe recent extensions to these programs, which 
will be the major emphasis of this paper are: 

(a) tc calculate the flow i*-. an inlet with a lead- 
in, enge slot, (b) to calculate the performance of 
suction and tangc.itiui blowing boundary layer con- 
trol concepts, and (c) to analyte the flow in aux- 
iliary inlet geoiik’tries at anglc-of -attack . Ihe 
mi'thod has also been applied to finding tne optimum 
velocity distribution in a subsonic diffuser. 

U'ometry Progr am 

A program called bllHCl is used for £-d and 
axisynmetnc geometries. For an axisymmetric inlet 
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case, the geometry Is represented by Its meridional 
profile which Is shown In Fig. 2(a). A i-o repre- 
sentation of the axisynmetric Inlet Is shown In 
Fig. 2(b). both the external and Internal aucts 
are extended far downstream so that accurate poten- 
tial flow solution can be obtained In the region of 
Interest. SCIACL breaks the profile Into segments 
with a control point on each senment which are used 
for potential flow calculations. The program also 
calculates Information such as curvature, wall 
angles, and flow area distrijutlon which are very 
useful In preliminary screening of proposed Inlet 
shapes. In addition to the surface points, SLl.tCt. 
generates oft body points (like flow measuring 
rakes) also shown In Fig. 2(a) at axial locations 
where the velocity profile or the streamlines are 
oesired. 

The J-d geometry program, applicable to Inlet 
geonetries like that shown In Fig. 2(c) and dis- 
cussed in Hef. 13, allows the user to Input a rela- 
tively small number of points to define the Inlet 
and centerbody. The routine then enriches the 
point number and redistributes the points fcr good 
potential flow analysis. The detailed description 
of this geometry package is given In Kef. b. borne 
examples of the geometries generated by this pro- 
gram are shewn in Figs. 2(b) and (c). 

Incompressible Po t ential Fl o w Basic bo lutlons 

The Oouglas Newman Programli^l^-lb Is used 
for calculating tne Incompressible potential flow 
field. The following basic solutions are obtained 
by tne above incompressible potential f Uiw program: 

1. btatic solution (V. • 0) 

2. Uniform axial flow solution 

3. go* angle of attack '.olution 

4. go’ angle of yaw solution (for 3-d geometry 
onlv) 

In general, to obtain the basic solutions, the sur- 
face IS replaced by a number of panels on which 
there Is a surface source (or sink) distribution of 
unknown strength. For 2-d and axisynmetric cases, 
the source density can be a constant, linear or 
parabolic. For the 3-d case, only a constant 
source density can be used. Ihe strength of source 
distribution varies over the surface In a manner 
such that at every control point the normal veloc- 
ity IS zero. However, the best static solution Is 
found to result from using a vorticity distribution 
on the cowl sui faces. 

Linear Lombination an d Corrections 

Ihe basic solutions obtained from the Incom- 
pressible potential flow calculation are combined 
linearly into a solution of interest having arbi- 
trary flo* conditions of free stream velocity, mass 
flow rate, and angle of attack.^^ In cases where 
a Xutta .ondition is required, the constants for 
linear combination are readjusted to satisfy the 
Kutta condition. Ihe linearly combined incompres- 
sible solution Is then corrected for compressibi I- 
ity.Tl* If the local velocity Is supersonic. It 
i« further corrected by the empirical supersonic 
correction formula. Ihe final potential flow 
solution can now be used as an input to the bound- 
ary layer program. 

Boundary Layer 

Ihe analysis of tne boundary layer uses a 
2-dimens1onal compressible boundary layer program. 


The complete documentation of the boundary layer 
program Is given In Ref. 6. The program calculates 
1 -.'•portant boundary layer parameters such as dis- 
placement thickness, momentum thickness, and skin 
friction coefficient, Cf. It also provides the 
boundary layer velocity profiles at any desired 
station. Tne location of transition from laminar 
to turbulent flow can either be predicted by the 
program or cai be specified by the user. Flow sep- 
aration Is defined to occur when the skin friction 
coefficient becomes zero. 

Recen. Extensions 


The discussion thus far has described th- ba- 
sic method of analysis. Now the discussion will 
focus cn describing the recent extensions which 
were motivated In part by the following thoughts. 

It Is les'rable to design a V/STOL Inlet as short 
and as thin as possible In order to reduce the 
weight, to reduce the friction drag at cruise, and 
Improve pilot visibility. However, when an Inlet 
Is too thin the peak velocity Is so high that the 
subsequent adverse pressure gradient causes the flow 
to separate at the lip resulting In a low pressure 
recovery and high distortion. Several ways to help 
control this possible separation are by an Inlet Up 
slot, the use of auxiliary Inlets, by suction O'* 
blowing boundary layer control, or by optimizing the 
surface pressure distribution. The analysis tech- 
niques to analyze these possibilities are considered 
next starting with the slotted Inlet. 

Slotted Inlet 

A slotted Inlet Is shown In Fig. 3. Two cases 
are considered, zero angle of attack and angle or 
attack. 

At zero angle of attack or at zero forward ve- 
locity for an axlsymtetrlc geometry the flow is ax- 
1 ymmetric and the Kutta condition (I.e., that the 
f.ow leaves the trailing edge of the slot, point 1 
on fig. 3, at a flow angle that bisects the slot 
trailing edge angle) Is applicable around the en- 
tire circumference of the inlet. Calculations were 
made fo'* stat'c :ond1tions, Vq • 0. Experimental 
data are Included tor comparison. The agreement Is 
quite good on the main Inlet cowl surface (points 6 
to 9). Ihe agreement between the theoretical and 
experimental surface velocities are good on the 
leading edge (points 4 and b) and the highlight 
ipoint 3). However, near the trailing edge (points 
1 and 2) the theory does not agree as closely with 
the experiment suggesting some modification to 
mathematical Kutta condition may be appropriate. 
Figure 3 shows that the peak velocities occurs at 
point 2 on the slat and point 7 on the main cowl, 
doth peaks are considerably lower than the peak ve- 
locity of the inlet without the slot,*'* also In- 
dicated on the figure. Inus the addition of the 
slot has unloaded the lip of the thin Inlet. 

Ihe case of slotted Inlet at an angle of at- 
tack IS more difficult, because the kutta condition 
can not correctly be Imposed at all circumferenllal 
positions simultaneously. 

Where the Kutta condition at an angle of at- 
tack IS required. It can be applied by adjusting 
the mass flow rate through the slot at one circum- 
ferential location until the Kutta condition Is 
satisfied at that circumferential location at the 
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trilling edge of the slat. Calculations were made 
using this approach ar.d the results are shown in 
Fig. 4. Experimental data are included for compar- 
ison. 

The Mach number was calculated at the circum- 
ferential ocstion of 270*. Agreement between 
theory and experiment is considered very good. 

Extending the method of analysis to include 
the Kutta condition results in the ability to ana- 
lyze a new class of inlet geometries, specifically 
those that employ leading-edge slats and slots. 

Auxiliary Inlet 

The method of analysis has also been extended 
to include auxiliary inlet geometries. Auxiliary 
inlets increase the total inlet flow area thereby 
reducing the amount of airflow that must be taken 
into the main inlet. It is another technique for 
preventing flow separation on cowl lips at static 
and low flight speed conditions. While an impor- 
tant application of auxiliary inlets is to super- 
sonic inlets at low speed, the application consid- 
ered here is to an auxiliary inlei on the top of a 
conventional subsonic inlet. The inlet and its 
panelling are shown in Fig. 5. A continuous N-line 
(ill the longitudinal direction) is required for the 
currant version of the 3-d potential flow program. 
When the N-line meets the auxiliary inlet opening, 
it is rerouted along the side wall of the auxiliary 
inlet and then proceeds back to the original N-line 
as shown in Figs. 5(a) and (b). Additional N-1 tnes 
are added to completely panel the inlet. This par- 
ticular example required 682 panels to describe the 
geometry. The 3-dimensional incompressible poten- 
tial flow code was then used to calculate the basic 
solutions at the center of each panel (control 
point) . 

With the ai .iliary inlet, the following tech- 
nique was found to yield the best static solution: 
Two inlet-duct systems are considered, one with a 
straight duct and one with a flared duct as shown 
in Figs. 5(a) and (b). The flared duct induces 
more flow through the main and auxiliary inlets. 

The difference between the velocities for the 
flared inlet and nonflared duct for a free stream 
uniform flow then provides the static solution. 

This procedure was adopted because the velocities 
in the region of an auxiliary inlet were unrealis- 
tically large when the vorticity distribution, 
noted earlier, was used for the static solution. 

The computer time for the basic solutions with 
the 682 panels is quite high, 19 minutes. However, 
the basic solutions are only computed once and are 
stored in the computer for later use in obtaining 
solutions of interest. Subsequent calculations 
using a linear combination method required only 
5 seconds of computer time. 


duced. For this case, the highest velocity occurs 
at the highlight. These sample cases indicate that 
this program can be used to calculate the surface 
velocities for nacelles employing auxiliary inlets 
and can pinpoint the problem areas. Figure 7 shows 
the effect of an auxiliary inlet on the peak veloc- 
Uy on the windward plane of the inlet. The peak 
velocity ratio is reduced from 2.8 to 2.5 at the 
highlight when the auxiliary inlet is opened. This 
is the desired result. 

Suction and Tangential Blowing 

Another recent extension to the basic methods 
is the analysis of suction and tangential blowing 
boundary layer control systems. Suction controls 
the boundary layer by removing that portion of it 
not having sufficient momentum to negotiate the 
subsequent adverse pressure gradient. Blowing con- 
trols the boundary layer by reenergizing it with a 
thin jet of high velocity air injected tangentially 
into the boundary layer. 

Some results from this analytical method are 
shown in Figs. 8 and 9 for suction and blowing, 
respectively. An axisymmetric inlet having a diam- 
eter of 0.508 m was analyzed at a free stream Mach 
number, M. - 0.12, throat Mach number, Mt - 0.4, 
and angle of attack, a • 60*. The skin friction 
coefficient distribution on the internal surface of 
the windwaro cowl is shown al.ng with boundary lay- 
er velocity profiles at several locations. Without 
boundary layer control, the solid line, the flow 
sep 'rates at S/L « U.81 where the skin friction 
coetiicient becomes zero. The boundary layer pro- 
file, just before separation, is quite weak com- 
pared to the one upstream at S/L - 0.48. Separa- 
tion is prevented when the boundary layer is con- 
trolled by suction (Fig. 8) - the dashed line. It 
was necessary to bleed off only 0.12 percent of the 
inlet lass flow to prevent separation as indicated 
by the nonzero skin friction coefficients. The 
static-to-total pressure ratio at suction location 
is 0.796. 

For the blowing boundary layer control (Fig. 

9) a blowing velocity ratio, jet velocity to 
boundary layer edge velocity, Vj/Vg - 1.75 was 
selected, for this case a blow mass flow of 0.4 
percent of the inlet mass flow was required to 
maintain attached flow. The reenergized boundary 
layer is clearly evident in the velocity profile 
just downstream of the blowing slot. 

Optimum Oiffuser Velocity Distribution 

Another application of the method is concerned 
with finding the optimum velocity distribution in a 
subsonic diffuser. This velocity distribution will 
result in the shortest no-boundary layer control 
inlet 'nd the lowest loss for the required amount 
of dif jsion. 


The velocity distribution along an N-line is 
illustrated in Fig. 6 for Vj;/V„ « 1.5 and 

1. a - 0*. 8 - 0* 

2. « - 30 , 8 - 0* 

The location is indicated by the letters along the 
N-line. At a • 0*, 8 ■ 0*, the peak velocity oc- 
curs at point 0 (Fig. 6) which is close to the cen- 
ter line on the downstream side of the auxiliary 
inlet, (lenerally speaking, the downstream surface 
of the auxiliary inlet is the high velocity area as 
mi^ht be expected. As angle of attack increases to 
30 , the peak velocities at 0 is substantially re- 


The mith^J of design of an optimum subsonic 
inlet is given in Refs. 20 and 21. Based on the 
design criteria given in those references, the 
boundary layer program was used to find the optimum 
diffuser velocity distribution. The generalized 
mathematical form of the velocity distribution is 
given by 


V * 
de 


<''m^x - ''de>* 


- 51 *’ 
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A typical case of Vn.jx ■ "’/^i ^de * 

shown In Fig. 10. The overall diffusion ratio 
Vjg Is the same for the three cases shown. The up- 
per part of the figure shows the surface velocity 
ratio (V/Vjp' as a function if surface distance 
s/Sppf. Velocity dlstributlins were calculated 
for three valu’S of the exponent b. 

A value of b - 0.613 produces the steepest 
Initial velocity gradient (largest Initial adverse 
pressure gradient). The Initial adverse pressure 
radlent Is so large that the flow separates on the 
Ip at the beginning of the diffusion process. A 
value of b > 1.006 produces a relatively more sev- 
ere adverse pressure gradient In the diffuser and 
the flow separates there. Somewhere between these 
two cases, there exists a velocity distribution 
such that at every location, the momentum of the 
boundary layer Is just able to overcome the adverse 
pressure gradient so that the flow remains attached 
throughout the diffuser. Ihls Ij called the opti- 
mum diffuser velocity distribution and Is achieved 
when b • 0.794. For comparison, Stratford's opti- 
mum velocity distributional is also presented in 
Fig. 10. The present optimum ve* ;i*. distribution 
Is more conservative at the liegi, i of the diffu- 
sion process than that of Siratfoiu. Since Strat- 
ford*' distribution is derived on the basis of zero 
skin .riction throughout the diffuser, it can be 
consi'lered as a limiting case. A design velocity 
distribution (besides having a safety margin 
against separation) should have a slightly more 
gradual start to the pressure rise (the decelera- 
tion of velocity) than that of Stratford. The 
present optimum velocity distribution can be a use- 
ful design approach. 

Concluding Re m arks 

An analysis method based on incompressible 
potential flow corrected for compressibility was 
described. Several sample calculations compared 
well with experimental data. The most recent ap- 
plications include inlet with a leading edge slot, 
an auxiliary inlet, and suction or blowing boundary 
layer control. An optimum diffuser velocity dis- 
tribution was also developed. This paper shows 
that the present methods can be a very powerful 
tool for the analysis of flow about and the design 
of V/STul inlets. 
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Figure 1. - Schematic diagram of computer programs. 
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Figures. - Auxiliiry Inlets. 




Figure 6. - The surface velocity of an auxiliary inlet 









